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ABSTRACT 

We have observed an unusually strong X-ray burst as a part of our regular 
eclipse timing observations of the low mass binary system EXO 0748— 676. The 
burst peak flux was 5.2 x 10~ 8 ergs cm _2 s _1 , approximately five times the normal 
peak X-ray burst flux observed from this source by RXTE. Spectral fits to the 
data strongly suggest that photospheric radius expansion occurred during the 
burst. In this Letter we examine the properties of this X-ray burst, which is 
the first example of a radius expansion burst from EXO 0748— 676 observed by 
RXTE. We find no evidence for coherent burst oscillations. Assuming that the 
peak burst luminosity is the Eddington luminosity for a 1.4 M Q neutron star 
we derive a distance to EXO 0748— 676 of 7.7 kpc for a helium-dominated burst 
photosphere and 5.9 kpc for a hydrogen-dominated burst photosphere. 

Subject headings: stars: distances — stars: individual (EXO 0748-676) — X-rays: 
binary — X-rays: bursts 
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1. Introduction 

The low mass X-ray binary (LMXB) system EXO 0748— 676 was originally discovered 
by EXOSAT during February 1985 (Parmar et al. 1985). Observations by EXOSAT de- 
tected periodic full eclipses and 25 Type I X-ray bursts, three of which were reported to 
show evidence of photospheric radius expansion (Gottwald et al. 1986). EXO 0748— 676 was 
subsequently observed by GINGA, ROSAT, ASCA, Beppo-SAX, XMM-Newton, and the 
Rossi X-Ray Timing Explorer (RXTE). EXO 0748-676 is a 3.82 hr orbital period system, 
normally classified as an "atoll" source, with a neutron star primary accreting matter from 
the Roche lobe-filling, low-mass main-sequence secondary star UY Vol. A possible observa- 
tion of gravitationally red-shifted lines in X-ray burst spectra has been reported by Cottam, 
Paerels, & Mendez (2002). Villarreal & Strohmayer (2004) have recently reported a possible 
spin period for the neutron star of 45 Hz. 

EXO 0748— 676 is one of the most heavily studied and important LMXB physics labo- 
ratories so constraining the source distance as accurately and rigorously as possible is very 
important. The spectral type of the secondary in EXO 0748— 676 has, so far, been elusive. 
An upper limit for the U— brightness of the secondary star of 21 st magnitude during X-ray 
quiescence can be estimated based on the absence of detectable emission on the SRC-J survey 
plate of the EXO 0748-676 region (Wade et al. 1985). Observations of EXO 0748-676 after 
its discovery by EXOSAT during its X-ray active state show an optically erratic variable of 
near 17 th magnitude (van Paradijs, van der Klis, & Pedersen 1988). However, van Paradijs 
et al. attributed the bulk of the optical emission to the neutron star accretion disk. The 
intrinsic faintness of the secondary and the system location outside the galactic bulge at 
(1, b) = (280.0°, —19.8°), makes an accurate reddening distance difficult to obtain (Schoembs 
& Zoeschinger 1990). Furthermore, its galactic latitude puts it well out of the galactic plane 
making a galactic differential rotation distance uncertain by an unknown amount (Cramp- 
ton et al. 1986). It is also too distant for the present generation of astrometry satellites 
to measure an accurate parallax. Yet the fact that EXO 0748— 676 is a bursting LMXB 
opens the possibility of obtaining the distance if an Eddington luminosity-limited X-ray 
burst can be observed and accurately analyzed. Gottwald et al. (1986) analyzed EXOSAT 
observations of three apparent radius expansion bursts but only assumed a distance of 10 
kpc to EXO 0748— 676 in order to estimate the Eddington luminosity. Jonker & Nelemans 
(2004) used the Gottwald et al. peak burst flux from the EXOSAT Medium Energy detector 
in order to estimate the distance of EXO 0748— 676 to be in the range 6.8-9.1 kpc. Until 
now, however, no X-ray burst detected by RXTE from EXO 0748— 676 showing evidence of 
photospheric radius expansion has been reported. 

Since 1996, near the beginning of the RXTE mission, we have maintained a program 
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of regular observations of EXO 0748— 676, of length 2-3 ks, in order to monitor the eclipses 
and precisely time the system orbital period (Hertz et al. 1997; Wolff et al. 2002). During 
our eclipse monitoring observations with the Proportional Counter Array (PCA) on RXTE 
we serendipitously capture numerous Type I X-ray bursts. A typical X-ray burst will rise 
within seconds to an apparent count rate of roughly 1000 counts s _1 PCU _1 in the 2-20 keV 
energy range (~ 1/3 Crab). However, during May 2004 we observed one X-ray burst that 
reached a count rate of nearly 5000 counts s _1 PCU _1 (~ 1.7 Crab) and showed evidence 
of photospheric radius expansion. In this paper we report a detailed analysis of this very 
luminous burst. Assuming the peak burst luminosity is equal to the Eddington luminosity 
for a 1.4 M Q neutron star, we derive a distance to the EXO 0748— 676 system of 7.7kpc 
for a helium-dominated burst photosphere, and 5.9 kpc for a hydrogen-dominated burst 
photosphere. 



2. RXTE Observations 

During RXTE observation designated by the ObsID 90059-03-01-00, we observed a 
strong X-ray burst that began at MJD (TT; not barycentered) = 53126.40358. Four PCUs 
(PCU0, PCU1, PCU2, PCU3) were active during this observation. The observation was 
done in GoodXenon mode which preserves full time and spectral information for each photon 
event. The burst peak flux saturated the GoodXenon data collection buffer twice during the 
burst, triggering the burstcatcher mode, which preserves only the event times binned at 122 
/is intervals and no spectral (PHA channel) or layer information. Time-dependent spectral 
analysis of the burst can be performed only in sections of the observation where there are 
no gaps in the GoodXenon data. However, we can construct a continuous light curve for the 
entire X-ray burst using Standard-1 data, which has no gaps. 

For the spectral analysis of this burst layer 1 events were extracted from the PCA FITS 
files for which unsaturated GoodXenon event data is available. We utilize version 5.3.1 of 
the HEASOFT FTOOLS package with XSPEC version 11.3.1 for the spectral fitting. This 
version of the FTOOLS generates response matrices with updated values of the default 
geometric areas for each PCU (Jahoda et al. 2005). These new effective areas improve the 
flux calibration of the Crab pulsar plus nebula so that the Crab flux is consistent with the 
standard value. This makes adjustments to the flux such as those utilized by Kuulkers et al. 
(2003) no longer necessary. Deadtime corrections are applied (always less than 7.6%) to the 
observed count rates and spectra. 

We fix the persistent X-ray emission for the source by extracting the spectrum of the 
200 seconds of data before the X-ray burst began and use this as a fixed background for the 
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burst spectral analysis. The Standard- 1 mode light curve is shown in Figure 1. Examination 
of Figure 1 shows that this pre-burst segment of data has relatively little flickering or dipping 
behavior that might signify locally intervening material in the EXO 0748—676 system and 
we therefore ignore this effect in our analysis. We fit the GoodXenon spectral data derived 
from segments of the burst profile to a model consisting of the XSPEC blackbody emission 
function multiplied by a photoelectric absorption function: wabsxbbodyrad. The interstellar 
hydrogen column density is frozen at 4.0 x 10 21 cm" 2 obtained by Thomas et al. (1997) 
from ASCA measurements. This allows the fitted parameters for the X-ray burst itself to 
vary independently of the persistent emission. For our spectral fits we include counts in 
the energy range 2.0—20 keV. We generate one response matrix for the entire X-ray burst 
observation which corresponds to approximately 100 seconds of time, and we assume zero 
systematic error in the fitting procedure. The resulting reduced x 2 values for the spectral 
fits are all in the range 0.5 — 2.0. This procedure is similar to the analysis method of Smale 
(2001) for an X-ray burst observed from X2127+119 in the globular cluster M15. As a test 
of our method, we analyzed the RXTE archival data for the same X-ray burst observed by 
Smale (2001) and found results similar to his, which confirms that our analysis procedure is 
reliable. 



3. Burst Oscillation Search 

The unusual strength of this burst makes it particularly attractive for burst oscillation 
searches. Villarreal & Strohmayer (2004) found evidence for a 45-Hz pulsation frequency in 
an incoherent sum of power spectra from 38 bursts from EXO 0748—676 but oscillations have 
never been detected in a single burst. We searched both the GoodXenon and burstcatcher 
data for coherent oscillations during this strong burst. 

The burstcatcher data were binned at 122 fj,s intervals for 10 s starting about 1 s before 
the beginning of the burst. These data are continuous, but have no energy resolution. The 
GoodXenon data have some gaps during the brightest part of the burst because the telemetry 
buffers filled. However, full energy resolution was available and we searched 100 seconds of 
data beginning at 0.5 s before the burst. We created three time series to search for pulsations, 
one from the burstcatcher data with 122 /is bins and two from the GoodXenon data with 
244 fj,s bins with energy selections of 2.5-6.0 keV and 6.0-60.0 keV. All three time series were 
searched for oscillations using the powspec FTOOL. A third-order polynomial was removed 
from each data window to reduce the red noise from the overall burst profile, and power 
spectra were Leahy normalized. Several combinations of FFT window length (from 2048 
to 32768 points) and numbers of power spectra averaged (from 1 to 10) were searched to 
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maximize sensitivity to oscillations that occur with various timescales. 

No oscillations were found during the X-ray burst. The precise upper limit on the 
oscillation amplitude depends on the assumed duration of the oscillations, but for averages 
of four 8192-point (2 second) intervals in the 6.0-60 keV time series, we would easily have 
detected sinusoidal pulsations whose RMS was at least 5% of the total RMS fluctuations. 
However, we would not have detected burst oscillations with an amplitude of only 3% such 
as those found by Villarreal & Strohmayer (2004). 



4. Burst Spectral Analysis 

In our discussion of the burst physics we adopt the notation of Ebisuzaki, Hanawa, 
& Sugimoto (1984), and define the observed "blackbody radius" i?bb with respect to the 
bolometric luminosity L and the observed color temperature T c using 

L = AnRl h aT*. (1) 

Solving for the blackbody radius yields 

/ T \ 1/2 

In the Ebisuzaki et al. (1984) development, the color temperature observed at infinity T c is the 
physical temperature of the radiating plasma at the "color radius" r c , where thermalization 
of the radiation occurs via either free-free absorption or Compton scattering. In this sense, 
the color radius can be interpreted as the standard photospheric radius. On the other hand, 
i?bb is n °t a true physical radius but rather a parameter that ensures that the observed 
spectral flux is correctly fitted in the bbodyrad model. Ebisuzaki et al. showed that i?bb is 
related to the color radius via r c = S Rbb, where S = (Stc/A) 1 ^ 2 is the "radius expansion 
factor" and r c is the electron scattering optical depth to the color radius (see also London, 
Taam, & Howard 1986; Smale 2001). Physically, the distinction between the two radii arises 
as a consequence of electron scattering. The presence of a relatively deep scattering layer 
above the color radius reduces the diffusion velocity of the photons, which in turn reduces 
the outgoing flux. In order to carry the required luminosity, r c must therefore exceed the 
value of i?bb computed using equation (2), since this expression completely ignores the effect 
of diffusion by implicitly assuming that the radiation is directly emitted from the "naked" 
stellar surface. Hence i?bb is significantly smaller than the true physical radius r c during the 
scattering-dominated, brightest part of the burst. This point is further discussed below. 
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The value of the Eddington luminosity LecM observed at infinity is given by the standard 
expression as modified by Lewin, van Paradijs, & Taam (1993) and Galloway et al. (2003), 



where M* is the neutron star mass (taken to be 1.4 M ) and k sc is the electron scattering cross 
section. The temperature-dependent factor in brackets in equation (3) expresses the high 
temperature modification of the scattering cross section. This set of general relativistically- 
correct relations allows us to determine the distance to EXO 0748—676. 

The X-ray burst takes ~ 2.5 s to rise to its peak count rate in the light curve shown 
in Figure 2. This rise takes the observed flux from a persistent source-plus-background 
count rate of 276 counts s -1 in the 2-20 keV energy range (~ 16 mCrab) to a burst peak 
rate over 18,000 counts s -1 . The EXO 0748— 676 orbital phase at the start of the burst is 
= 0.034220(4), roughly 230 s after eclipse egress. After this the flux spends ~ 1.7 s within 
90% of the peak count rate. The flux then falls rapidly for ~ 2 s down to a plateau at a level 
~ 30% of the peak rate and then starts a rapid exponential decay with an e-folding time of 
14 s. The principal parameters we derive for this burst are given in Table 1. 

In Figure 2 the results of our spectral fitting are presented for the helium-rich photo- 
sphere case. The hardness ratio we plot is the ratio of the flux in the 6-20 keV energy band to 
the flux in the 3-6 keV energy band. When the observed flux reaches its peak, the apparent 
blackbody radius i?bb also reaches its maximum value as seen in Figure 2. As the observed 
blackbody radius initially expands to ~ 20 km, the observed color temperature decreases to 
a local minimum at the moment the radius achieves its largest value. In a similar fashion, the 
hardness ratio falls reaching a local minimum as the observed flux is reaching its peak. The 
hardness ratio then dramatically increases again reaching its maximum value ~2 seconds 
after the flux has begun to decay. The decrease in both the temperature and the hardness 
ratio while the flux is near its peak is indicative of the cooling that occurs during the photo- 
spheric expansion (Strohmayer & Bildsten 2003). After the flux reaches its peak and begins 
to decay, the radius drops rapidly to an apparent minimum of ~4.6 km and then recovers 
and asymptotically approaches ~13 km. As i?bb decreases after the peak, the temperature 
quickly rises to its largest values just as the flux is also decreasing. When the radius reaches 
its minimum value, just after the peak of the burst, the temperature and the hardness ratio 
each reach their respective maxima. After this the burst begins its exponential decay until 
it fades away completely. 

We derive an unabsorbed burst fluence of 3.60 x 10~ 7 ergs cm -2 by multiplying the 
unabsorbed flux in the energy range 2.0-20.0 keV by a bolometric correction and then inte- 
grating over the duration of the burst rise and decay. The bolometric correction is computed 
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as a function of time by dividing the total blackbody flux by the blackbody flux in the 
2.0-20.0 keV energy range based on the sequence of observed color temperatures. Assuming 
that the expanding envelope is helium-dominated, and that the peak flux of 5.16 x 10~ 8 
ergs cm -2 s _1 corresponds to the Eddington luminosity of 3.64 x 10 38 ergs s _1 for a 1.4 M 
neutron star, we derive a distance of 7.67 kpc to EXO 0748— 676. On the other hand, if the 
burst photosphere has solar abundances, then the peak flux corresponds to an Eddington 
luminosity of 2.12 x 10 38 ergs s _1 for a 1.4 M neutron star, and a distance of 5.85 kpc to 
EXO 0748—676. The associated values of L E dd for hydrogen- and helium-rich photospheres 
are consistent with the peak luminosities obtained by Kuulkers et al. (2003) in their anal- 
ysis of photospheric radius expansion bursts in globular cluster sources, which have known 
distances. We note that the peak flux we obtain is ~ 26% higher than the apparent peak 
flux obtained by Gottwald et al. (1986) for three photospheric radius expansion bursts from 
EXO 0748-676 observed by EXOSAT. Galloway et al. (2003) found that in photospheric 
radius expansion bursts from the LMXB 4U 1728—34 the peak fluxes were not constant but 
instead exhibited variations with a RMS deviation of 9.4% and a total variation of 46%. 
Galloway et al. concluded that this variation might be caused by variable obscuration from 
a warped accretion disk around the neutron star. If the accretion disk surrounding the neu- 
tron star in EXO 0748—676 is also warped, then this could cause the same type of variations 
in peak fluxes as observed in 4U 1728—34, which could explain why our peak flux is higher 
than that of Gottwald et al. (1986). If our higher flux is a result of less obscuration by a 
warped accretion disk then the distance we derive will be more accurate than a distance 
based on the lower EXOSAT flux because our peak flux is more indicative of the unobscured 
Eddington flux for EXO 0748-676. 

The EXO 0748— 676 system is located at galactic latitude —19.8° and galactic longitude 
280.0°. Hence if we assume that the galactic center is 8.5 kpc from the solar system, then 
EXO 0748—676 is ~ 10.5 kpc from the galactic center and ~ 2.6 kpc below the galactic plane. 
The ratio of the persistent X-ray flux from the pre-burst spectrum to the Eddington flux at 
the peak of the burst is 7 = 0.0094 (see Table 1), which is consistent with the assumption 
that EXO 0748— 676 is an atoll LMXB source. The distance implied by the helium-rich 
photosphere case combined with the observed persistent X-ray flux yields a mass accretion 
rate of ~ 2.9 x 10~ 10 M yr" 1 , which is in the accretion rate range expected to yield mixed 
hydrogen and helium burning triggered by thermally unstable helium ignition (Strohmayer 
& Bildsten 2003). On the other hand, if the burst photosphere is hydrogen-rich then the 
closer distance implies a lower persistent M, pushing M into the range that results in burst 
ignition via thermally unstable hydrogen burning (Strohmayer & Bildsten 2003). Thus, 
conclusions about the characteristics of the burning for this X-ray burst depend on the 
assumed composition of the burst photosphere. 
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A number of possible sources of error exist for our estimated distance to EXO 0748—676. 
First, if the assumed mass of the neutron star (M* = 1.4 M Q ) is incorrect then the Eddington 
luminosity has been incorrectly calculated and the distance must be adjusted accordingly. 
For a 1.6 (1.2) M Q neutron star we derive a distance of 8.2 (7.1) kpc from the same spectral 
fitting data for a helium-dominated burst. A second source of error enters through our 
estimate of the burst bolometric flux, either due to errors in the PCA calibration or errors 
in the bolometric corrections derived from the fitted temperatures. If the flux is 20% too 
high (low) we estimate that the source distance would change to 7.0 (8.6) kpc, again for a 
helium-dominated burst. The error estimates for the hydrogen-rich burst case are similar. 

The increase in R hh during the early phase of the burst clearly suggests that the pho- 
tospheric radius of the radiating atmosphere expands. However, the subsequent minimum 
value of i?bb ~ 4.6 km just after the burst peak is only about one third of the apparent 
radius that a distant observer would associate with a 1.4 M neutron star with a proper 
radius of 10 km, which is about 13 km when the gravitational redshift is included (Shapiro 
& Teukolsky 1983; Lattimer & Prakash 2001). If one expects a proper neutron star radius 
of ~ 10 km, then that in turn suggests that the bbodyrad model in XSPEC does not include 
all of the physics involved in producing a radius expansion X-ray burst light curve. Without 
a clear understanding of the physics involved in the light curve production after the peak 
of the burst, the derived radii cannot be reliably used to constrain either the neutron star 
radius or the equation of state of the nuclear matter (see the discussion of this point in 
Lewin et al. 1993). During the early, high-temperature portion of the burst around the time 
of peak flux, the thermalization of the emitted radiation is almost certainly dominated by 
Comptonization. In this situation, the scattering depth to the color radius, r c , is related 
to the observed color temperature, T c , via r c = (m e c 2 /AkT c ) 1 / 2 (London et al. 1984, 1986; 
Smale 2001). Based on the observed color temperature at the peak of the burst, we find 
that r c ~ 5, and therefore S ~ 2, where S is the radius expansion factor defined above. This 
implies that the color radius r c is roughly double the blackbody radius just after the burst 
peak, when i?bb achieves its minimum value. However, even a factor of two increase in i?bb 
is not sufficient to bring the color radius r c into accord with the apparent radius of 13 km. It 
should be noted that the canonical value of 10 km for the proper radius of the neutron star 
is based on a particular equation of state for the nuclear matter, and some alternative equa- 
tions may allow slightly more compact neutron stars (Heiselberg & Pandharpande 2000), 
with proper radii perhaps as small as 9 km. However, this adjustment does not appear to 
be sufficient to completely remove the discrepancy cited above. 

With regard to the post-peak-flux drop in R hh it is possible that non-thermal spectral 
distortions could be created by ionized material located along the line of sight. Such ma- 
terial could rise from the accretion disk due to strong heating and photoionization driven 
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by the peak X-ray flux. Preferential absorption of the soft radiation by the ionized gas 
would decrease the total flux, while also causing the apparent color temperature to rise, in 
agreement with the observations (see Figure 2). These two effects will cause a reduction 
in the value of i?bb computed using equation (2). In this sense, the X-ray burst will have 
created its own "weather" that it must shine through if the photons are to reach the PCA. 
There is some evidence for the existence of ionized gas in the vicinity of the binary system 
4U 1820—303 during X-ray bursts although no detailed model for this phenomenon has ap- 
peared (Strohmayer & Brown 2002). Presumably, after the early phase of strong irradiation, 
the ionized material drops back down into the disk and therefore the differential absorption 
effect disappears. However, it is not clear if the dynamical and thermal timescales in the gas 
support this interpretation. 

5. Conclusions 

We have analyzed an unusually strong X-ray burst from the LMXB EXO 0748—676 and 
found that it shows evidence for photospheric radius expansion. No oscillations were detected 
during the X-ray burst despite an exhaustive search. When we equate the peak flux from 
the X-ray burst with the expected Eddington flux from a 1.4 M neutron star we derive a 
distance to the EXO 0748— 676 system of 7.7 kpc for a helium-dominated burst photosphere 
and 5.9 kpc for a hydrogen-dominated burst photosphere. While we have argued that our 
distance determination to EXO 0748— 676 is likely to be more accurate than previous deter- 
minations the ambiguity in the distance determinations between hydrogen- and helium-rich 
photospheric abundances in the real Eddington luminosity remains. Thus, we must quote 
a range of distances accounting for both abundance cases. Only when the observational 
characterizations of burst properties are able to fix the abundances of the bursting neutron 
star photosphere can this ambiguity in distance determinations to bursting LMXB systems 
be overcome. The X-ray burst flux light curve strongly suggests that we do not completely 
understand the physics of the spectral formation process in an X-ray burst just after max- 
imum brightness, when the luminosity is beginning to recede from the Eddington limit. If 
the derived distance is correct, then the small values obtained for the fitted blackbody radius 
i?bb after the peak of the burst are difficult to reconcile with the 10 km proper radius of the 
star, even when adjustments are made based on the high electron scattering optical depth. 
Possible alternative explanations include the presence of intervening material that absorbs 
the soft flux after the peak, or perhaps the occurrence of non-spherical mass motions. 

We are pleased to acknowledge discussions with Drs. Craig Markwardt, Jean Swank, 
Deepto Chakrabarty, Duncan Galloway, Jeroen Homan, and Lev Titarchuk. This research 
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TIME (s) 

Fig. I. — The Standard-1 light curve of the X-ray burst from the four active PCUs. The 
dashed line shows the approximate background count level in the PCA during the burst. 
The pre-burst light curve shows only a small level of variation consistent with little or no 
dipping behavior as noted in the text. 
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Fig. 2. — Results for the bolometric flux, observed color temperature, and blackbody radius 
from the spectral fit during the X-ray burst assuming a distance of 7.7 kpc. Note the gaps 
in the data points beginning at 1.6 s and 3.1 s. These correspond to the saturation of the 
GoodXenon event mode. The peak bolometric flux is F pcak = 5.16 x 10~ 8 ergs cur 1 s _1 . 
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Table 1. EXO 0748-676 Burst and Derived Parameters 



Parameter 



Value 



Burst Fluence (= E b ) 
Peak Bolometric Flux (= F pea k) 
Persistent Flux (= F pcr ) 

T [ -^per/ (-^peak -^per)] 

a [= F pcr AT/E b ] a 
r [= E b /F peak ] 



3.60 x 10 7 ergs cm 2 



5.16 x 10" 



ergs cm 2 s 1 

-2 -l 



4.81 x 10 10 ergs cm ' s 
0.0094 
17 
7.0 



Derived Parameters Assuming Hydrogen Burst Photosphere 



EXO 0748-676 Distance From Sun b ' c 
EXO 0748-676 Distance To Galactic Center b 
EXO 0748-676 Distance Below Galactic Plane 
Burst Energy Output 
Persistent Accretion Rate M per 



5.9 ±0.9 kpc 
9.5 kpc 
2.0 kpc 
1.48 x 10 39 ergs 
1.68 x 10- 10 M /year 



Derived Parameters Assuming Helium Burst Photosphere 



EXO 0748-676 Distance From Sun b ' c 
EXO 0748-676 Distance To Galactic Center b 
EXO 0748-676 Distance Below Galactic Plane 
Burst Energy Output 
Persistent Accretion Rate M per 



7.7 ±0.9 kpc 
10.5 kpc 
2.6 kpc 
2.54 x 10 39 ergs 
2.89 x 10~ 10 M /year 



a Assuming a burst recurrence time of AT = 0.285 hours (Galloway et al. 2005) 
b Assuming a solar system distance to the galactic center of 8.5 kpc. 
c Assuming L peak = L Edd . 



